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I. INTRODUCTION 



The phase transition from partonic degrees of freedom (quarks and gluons) in ultra- 
relativistic nuclear coUisions to hadronic degree of freedom is a central focus of recent ex- 
periments at RHIC [1-6]. One of the interesting and unexpected discoveries [7,8] at RHIC 
is the " baryon anomaly" [9] , observed as a large enhancement of the baryon to meson ratio 
and as a large difference between the nuclear modification factor Raa{pt) between total 
charged and 7r° at moderate 2 < px < 5 GeV/c. There are two main effects that contribute 
to this anomaly. One is the predicted jet quenching [9,10] that strongly suppresses the pion 
yield above pr > 2 GeV/c. This effect causes an apparent enhancement since the pion 
denominator is reduced. The second effect is a genuine enhancement of baryon transverse 
momentum spectra. Several effects can contribute to such enhancements of baryon yields at 
moderate pr- Radial hydrodynamical flow has been observed at all energies including RHIC 
[11]. However at high px local equilibrium must fail. At intermediate pt a nonequilibrium 
remnant of hydrodynamic flow may arise from multi quark recombination [6,12]. 

In this work we continue to explore the possibility that a more novel and unconventional 
source of baryon production may be at least partially responsible for the baryon anomaly. 
We study whether baryon junction loops, as proposed in [13] to explain (anti)hyperon pro- 
duction at lower (SPS) energies, could help explain the RHIC data. The idea that nonper- 
turbative three color flux junctions could play an important role in baryon and anti-baryon 
production at high energies was proposed long ago by Rossi and Veneziano [14,15] on the 
basis of dual regge theory. This idea was extended and applied by Kharzeev [16] to nuclear 
collisions. The flrst full A-l-A event Monte Carlo implementation of this mechanism was 
HIJING/B generator by Vance et al. [17]. The addition of baryon junction loop mechanism 
led to HIJING/BBl.lO version of this model [13]. Unlike conventional diquark fragmenta- 
tion models, a baryon junction allows the diquark to split with the three independent flux 
lines tied together at a junction. For an alternate possible formulation of baryon junction 
dynamics see [18-21]. 
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In this paper we introduce a new version (v2.0) of HIJING/BB that differs from HI- 
JING/BBl.lO [13] in its implementation of hypothesized junction anti-junction loops that 
may be responsible for novel baryon pair production in nucleus nucleus (A+B) collisions. 

A large data base on meson and baryon spectra [7,8,22-37] are now available from RHIC 
experiments. The HIJING event generator [38] was developed to extrapolate hadron-hadron 
multiparticle soft plus hard phenomenology as encoded in the LUND JETSET/PYTHIA 
model [39] to nuclear collisions. One important feature of HIJING is that it can account for 
the pion quenching component of the baryon anomaly. However, the LUND JETSET di- 
quark string fragmentation mechanism used in HIJING vl.37 [38] completely fails to describe 
the baryon spectra observed in A+A coUisions at all energies [40-42]. HIJING/BBl.lO was 
developed to address this failure at SPS energies. It was however also found to be also 
inadequate, as we review below, with respect to baryon observables. For a recent discussion 
comparing HIJING1.37 and HIJING/BBl.lO predictions for global observables at RHIC see 
[42]. 

Clarifying the physical origin of the (anti) baryon dynamics at RHIC is important given 
the variety of hadronization mechanisms proposed in hydrodynamic models [43], multi- 
quark coalescence [12], thermal models [44], and the novel baryon junction dynamics [9,45]. 
The valence baryon number migration over the large rapidity window —5 < y < 5 at 
RHIC provides another stringent tests of the baryon dynamics. Recent RHIC data show 
at midrapidity a sizeable finite net-proton (p — p) number in the final state [7,8,22-37]. 
Moreover, the antibaryon-to-baryon ratios are not equal to one, providing further evidence 
for non-transparency of high energy nuclei. This significant baryon number transport over 
more than 5 units in rapidity has inspired other approaches as well [46-50] . The net baryon 
density at midrapidity also has an impact on the hadro-chemical equilibration affecting 
hadronic yields [51]. 

Another possible source of novel baryon/hyperon physics are strong color electric fields 
(SCF). This is modeled as an increase of the effective string tension that controls the qq 
and qqqq pair creation rates [52,53]. Molecular dynamics model [54-57] have been used 
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to study the effects of color ropes as an effective description of the non-perturbative, soft 
gluonic part of QCD [58-60]. SCF increases the partons "intrinsic transverse momentum" 
{kr) and decreases the protons and antiprotons yields. The empirical value of the Regge 
slope for baryon is a ' ~ 1 GeV~'^ that yields a string tension k, ( related to the Regge slope, 
K = l/2TTa'), of approximately 1 GeV/fm. It has been suggested that the magnitude of 
a typical field strength at RHIC energies might reach 5-12 GeV/fm [61]. However, in this 
work we will not consider further the consequences of SCF effects but concentrate only on 
the junction loops. 

In the following sections we compare numerical results of HIJING/BB v2.0 to transverse 
momentum spectra, rapidity densities (dN/dy) of protons (p), anti- protons (p), and net- 
protons (p-p) for central Au-|-Au collisions at y^sj^=200 GeV as well as their centrahty 
dependence. The characteristic stopping observables: average rapidity loss, energy loss of 
net-baryons per participant nucleon and transverse energy per net-baryons are presented. 
The "anomalous" baryon-mesons composition at moderate pt observed in particle ratios 
p/p, p/tt"*", p/tt~ and in species dependent nuclear modification factors, is also discussed. 

In Section II we briefly recall the ingredients of the HIJING [38] , HIJING/BB (vl.lO) 
[17,13] approaches and point out the extensions incorporated in the new version HIJING/BB 
v2.0. In Section III numerical results are discussed in detail in comparison with available 
data. Summary and Conclusions are presented in Section IV. 

II. HIJING/BB V2.0. 

HIJING is a model that provides a theoretical framework to extrapolate elementary 
proton-proton multiparticle phenomena to complex nuclear collisions as well as to explore 
possible new physics such as energy loss and gluon shadowing [49]. Three versions of this 
model will be compared in this paper: HIJING vl.37, HIJING/BB vl.lO and HIJING/BB 
v2.0. Detailed descriptions of the HIJING vl.37 and HIJING/BB vl.lO models can be found 
in Refs. [13,17,38,42]. 
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In HIJING1.37 [38] the soft beam jet fragmentation is modeled by diquark-quark strings 
as in [39] with gluon kinks induced by soft gluon radiations. The mini-jets physics is com- 
puted via an eikonal multiple collision framework using pQCD PYTHIAS. 3 to compute the 
initial and final state radiation and hard scattering rates. The cross section for hard parton 
scatterings is enhanced by a factor K—2 in order to simulate high order corrections. HI- 
JING extends PYTHIA to include a number of new nuclear effects. Besides the Glauber 
nuclear eikonal extension, shadowing of nuclear parton distributions is modeled. In addition 
dynamical energy loss of the (mini) jets is taken into account through an effective dE/dx. 

In HIJING/BB vl.lO [13] the baryon junction mechanism was introduced as an extension 
of HIJING/B [17] in order to try to account for the observed longitudinal distributions 
of baryons(B) and anti-baryons(B) in proton nucleus (p+A) and nucleus-nucleus (A-l-A) 
collisions at the SPS energies. However, as implemented in HIJING/BBvl.lO the junction 
loops still fail to account for the observed transverse slopes of anti-baryons at moderate pr 
as shown in [42]. This motivated us to try to reformulate the junction loop implementation 
in the present HIJING/BB (v2.0). 

The J J loop algorithm of vl.lO has been replaced by a simple enhancement of the intrinsic 
(anti)diquark pt kicks in any string that has been selected to contain one or more loops. 
Multiple hard and soft interactions proceed as in HIJING1.37. Before fragmentation however 
via JETSET we compute the probability that a junction loop occurs in the string. A picture 
of a juction loop is as follows: a color flux line splits at some intermediate point into two 
flux line at one junction and then the flux line fuse back into one at a second anti-junction 
somewhere further along the original flux line. The distance in rapidity between these 
points is chosen via a Regge distribution as described bellow. For single inclusive baryon 
observables this distribution does not need to be specifled. 

The probability of such a loop is assumed to increase with the number of binary in- 
teractions, Uhits that the incident baryon suffers in passing through the oncoming nucleus. 
This number depends on the relative and absolute impact parameters and is computed in 
HIJING using the eikonal path through a diffuse nuclear density. 
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We assume as in [13] that out of the non single diffr active NN interaction cross section, 
Cin ~ CTnsd, & fraction fjj = (Jjj/{ain — (Tnsd) of the events excite a junction loop. The 
probability after nhus that the incident baryon has a J J loop is: 

Pjj^i-ii-fjjY--^ (1) 

We take ctj j= 17 mb, (Tsdf ~ 10 mb and the total inelastic nucleon nucleon cross section 
Gin ~ 42 mb at RHIC energies. These cross sections imply that a junction loop occurs in pp 
collisions at RHIC energy with a rather high probabihty 17/32 fa 0.5 and rapidly approaches 
1 in AA. In p + S where nhu ~ 2 there is an 80% probability that a junction loop occurs 
in this scheme. Thus the effects of loops is taken here to have a very rapid onset and 
essential all participant baryons are excited with J J loops in AA at RHIC. We investigate 
the sensitivity of the results to the value of parameter ajj and found no significant variation 
on pseudo-rapidity distributions of charged particles, for 15 mb < Ujj < 25 mb for Au+Au. 
Light ion reactions like p-\- S and p-\- Ar would have more sensitivity to (Tjj. 

The production of a baryon and antibaryon from a J J loop is simulated via an enhance- 
ment of the diquark pt kick parameter a^q = PARJ{21) of JETSET7.3. The default value 
is (Tqq — 0.36 GeV/c in ordinary string fragmentation. However in events where the string 
has a junctions loop we can expect a significantly higher pt kick [13]. We therefore propose 
a very simple algorithm whereby the J J is modeled by enhancing Uqq by a factor F^^ which 
we fit to best reproduce the observed pr spectrum of the baryons. This implementation of 
the J J model marks a radical departure from that implemented in HIJING/BBvl.lO. 

While the above model allows the baryon antibaryon pairs to acquire much high trans- 
verse momentum in accord with observation, the absolute production rate also depends on 
the diquark/quark suppression factor PARJ{1). The JETSET default for ordinary (funda- 
mental flux) strings has PARJ{1) — 0.1. The reduced number of protons and anti protons 
observed at RHIC relative to HIJING1.37, will be shown below to be accounted for, if 
PARJ(l) is reduced to 0.07 in J J loops. 

In summary two paramters, P ARJ {1), PAR J {21), are used in version 2.0 to simulate 
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the dynamical consequences of the hypothesized J J loop production m A -\- B reactions. 
The factor Fp^ modifying the default 0.36 GeV value of PARJ{21) may depend on beam 
energy, atomic mass number (A), centrality (impact parameter). However, we will show 
that a surprising good description of a variety of observables is obtained with a constant 
value, Fpy=3. The sensitivity of the theoretical predictions to this parameter is discussed 
in section III. 

Finally, we remark that correlations studies in p+p and p(d)+Au collisions at RHIC 
energies could eventually help us to obtain more precise values of J J loop parameters (mainly: 
(Tjj, Regge intercept q;j(0) and Fp^). The contribution to the double differential inclusive 
cross section for the inclusive production of a S and B in NN collisions due to J J exchange 
is [13,16]: 

a PbCI Pb 

where Cbb -is an unknown function of the transverse momentum and Mq + P + B 
(junction+Pomeron+baryon) couplings [13,17]. The predicted rapidity correlation length 
(1 — aj{0))~^ depends upon the value of the intercept aj{0). To test for Mq component 
aj{0) ~ 0.5 requires the measurement of rapidity correlations on a scale \yB — ys| ~ 2. In 
contrast, infinite range rapidity correlations are suggested if aj{0) ~ 1.0 [20]. It is thus 
important to look for rapidity correlations at RHIC energies where very high statistics data 
are now available, inp + p^B + B + X, or p{d) + A^B + B + X. 



III. NUMERICAL RESULTS 



A. Transverse Momentum Spectrum 



The nuclear modification factor (Raa) is defined as the ratio of the hadron yield in 
central Au+Au coUisions to that in p+p reactions scaled by the number of binary coUisions 

{Ncoii): 



7 



Raa{p±) 



(fNAA/dydp^ 



(3) 



< Ncoii > (PNpp/dydp± 



where, < Ncoii > is the average number of binary coUisions of the event sample calculated 
from the nuclear overlap integral (Taa) and the inelastic nucleon-nucleon cross section; 



In Fig. 1 the measured [62] nuclear modification factor {Raa) for charged hadrons in 
central (0-10%) Au+Au coUisons at 200 GeV are compared to the predictions of HIJING 
vl.37 and HIJING/BB v2.0 models. The data shows the strong jet quenching effect that 
suppresses the hadrons yield by a factor of ~ 5 for the highest px bins resulting in an observed 
maximum in Raa at pt ~ 2 GeV/c. Note that both HIJING vl.37 and HIJING/BB vl.lO 
fail to reproduce the "baryon bump" at moderate px seen in the Raa factor and also fail 
to account for the large transverse slopes of baryons and anti-baryons (see ref. [42] ). The 
Lund string fragmentation mechanism of hadronization in HIJING vl.37 leads to a rather 
slow increase of the nuclear modification factor Raa to unity at high px, not observed in the 
data (see part a, results without quenching and shadowing effects (nqs)). The addition of 
jet quenching and shadowing effects (yqs) in HIJING vl.37 still fails to describe the data. 

In contrast, HIJING/BB v2.0 with shadowing and jet quenching (yqs) effects with default 
energy loss parameter dE/dx=l GeV/fm (for quark jet) describes well the data over the full 
Pt range. Some of the observed discrepancies could be attributed to strange and multistrange 
hyperons which are underestimated in our calculations because we do not consider SCF 
effects here. 

Figure 2 presents a comparison of the experimental transverse mass distributions [22] 
of positive (left) and negative (right) particles with the predictions of HIJING/BB v2.0 
(upper panel) and HIJING vl.37 (lower panel). The data shows a mass dependence in 
the shape of the spectra. The protons (p) and anti-protons (p) spectra have a shoulder- 
arm shape at low pr characteristic of a radial flow. The pion spectra are well described 
by both models. Introducing corrected JJ loops algorithm in HIJING/BB v2.0 result in a 
significant improvement in the description of the protons and anti-protons in the scenario 



< Ncoii >= cr, 



.inel 
nn 



< Taa >. 
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with shadowing and jet quenching. However, only a quahtative description is obtained for 
low rriT — mo spectra due to the presence of radial flow, not included in the model. A similar 
conclusion can be drawn from the predictions of the models for mean transverse momenta. 

B. Stopping Observables 

Rapidity distributions of participants (net) baryons are very sensitive to the dynamical 
and statistical properties of nucleus-nucleus collisions. The RHIC net proton distribution is 
both quahtatively and quantitatively different from those at lower AGS and SPS energies 
[28]. Recent results for net-proton in central (0-5%) Au-|-Au interactions at total nucleon 
nucleon centre of mass energy y^Siviv=200 GeV show an unexpectedly large rapidity density 
at midrapidity [23,28]. 

Figure 3 presents a comparison of the rapidity distributions of protons (Fig. 3a), anti- 
protons (Fig. 3b) and net-protons (p-p) (Fig. 3c) and their ratio (p/p) (Fig. 3d) obtained 
for Au-|-Au interactions at total cm energy -v/s=200A GeV with the model predictions of 
HIJING/BB v2.0 and RQMD v2.4 [53]. Corrections for feed-down contributions have been 
applied to the data. We discuss here a comparison with RQMD v2.4 results in order to 
investigate if hadronic rescattering only and SCF effects as implemented in RQMD, could 
explain the new data. The new version of HIJING/BB reproduces very well the experimental 
yield at mid-rapidity for both p and p as well as their ratio and the net protons yield (p-p). 
This agreement is improved if shadowing and jet quenching are included. In contrast RQMD 
v2.4 does not reproduce the shape of the proton rapidity distribution near midrapidity and 
strongly underpredicts the anti-proton yield. 

In addition, centrality dependence of proton and anti-proton yields at mid-rapidity have 
been also analysed and the results are shown in Fig. 4. HIJING/BB vl.O overpredicts the 
data [23] except for very peripheral collisions. In contrast, HIJING/BB v2.0 reproduces very 
well the experimental yield at all centralities. 

One of the main feature of the data is the observed increase of net-proton up to three 
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units of rapidity away from midrapidity (Fig. 3c). This central valley could be used as an 
indicator for partonic processes [63] , [64] . Microscopically, the baryon number transport over 
4-5 units of rapidity to the equilibrated midrapidity region is not only due to hard processes 
acting on single valence (di) quark that are described by perturbative QCD, since this yields 
insufficient stopping [40], [41]. Instead, additional processes such as the nonperturbative 
junction mechanism as implemented in HIJING/BB v2.0 are able to reproduce the observed 
distribution. Such mechanism may lead to substantial stopping even at LHC energies. 

The net-baryon (B-B) distribution retains information about the energy loss and allows 
the degree of nuclear stopping to be determined. Experimentaly, to obtain the net-baryons, 
the number of net-neutrons and net-hyperons have to be estimated. In addition, the data 
need to be extrapolated to full rapidity space introducing other systematic errors. In con- 
trast, in models we can calculate directly specific stopping observable as the average rapidity 
loss and the energy loss per participant nucleon as defined in ref. [28] . The average rapidity 
loss is defined as < 5y >— Up— < y > where Up is the rapidity of the incoming projectile 
and < y > is the mean net-baryon rapidity after the collision: 



where Npart is the number of participating nucleons in the colhsion. 

The total energy Etot per net-baryon after the collisions can be derived using the relation: 



where < rriT >—< yjpr + > is the average transverse mass. The energy loss per partic- 
ipant pair could be estimated as AE — (100.0 — Efot) GeV. The predictions of the models 
for these quantities are presented in Fig. 5. 

Fig. 5a shows the net-baryons model predictions in comparison with BRAHMS data [28] 
obtained from the measured net-proton {N(^b_^^ fa 2A^(p_p)). At RHIC energies a broad 
minimum has developed at mid-rapidity for net-baryon spanning few units of rapidity indi- 
cating that collisions are quite transparent. The average net-baryons rapidity loss deduced 




(4) 




(5) 
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by BRAHMS < Sy >= 2.0 ± 0.2 [28] is well reproduced by HIJING/BB v2.0. In contrast, 
the experimental value for the total energy per net-baryon after the collisions, E^^-f = 28 ± 6 
GeV [28] , are far from our theoretical predictions of both HIJING models, Etot ~ 40 GeV. 

The RQMD model [53] with rescattering and SCF effects strongly overpredicts the mid- 
rapidity net-baryon distributions and the average rapidity loss. Partly as a consequence Efot 
per baryon after the collisions is underpredicted. The extrapolation to high rapidity used in 
[28] to obtain their values oi < 5y > and Etot may explain part of the observed discrepencies. 
A precise measurement of transverse energy per baryon (as shown in Fig. 5c) could help 
also in the study concerning the origin of rescattering (which could influence the dynamics 
of the reaction at hadronic or partonic stage). However, futher analysis and baryon rapidity 
distribution measurements at large rapidity are needed in order to draw a final conclusion 
and to use these observables as a signature for partonic processes and for quark-guon plasma 
(QGP) formation. 

C. Particle Ratios versus pT 

In Figure 6 the p/vr" ratio are shown as function of transverse momentum (pt) for central 
(0-10 %) (upper part) and peripheral (60-90 %) coUisions (lower part). The measured ratios 
increase at low pr and saturate at different values which strongly increase from peripheral 
to central collisions. For central collisions at moderate pt (2 < Pt < 5 GeV/c ), the yields 
of anti-protons (p) are comparable to that of pions. In hard scattering processes described 
by pQCD and implemented in HIJING, the ratio p/tt" are determined by the fragmentation 
of energetic partons, independent of the initial colliding system. Therefore within the errors 
those ratios are well described for peripheral coUisions by both version of the models. As 
expected in peripheral collisions there is little sensitivity to the new ingredients implemented 
in HIJING/BB v2.0. In contrast the clear increase in the p/tt" ratio at moderate pt from 
peripheral to central collisions, is seen to be sensitive to the new physics as implemented 
in HIJING/BB v2.0. HIJING vl.37 strongly underpredicts the observed ratios at moderate 
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Pt, an effect that is corrected in HIJING/BB v2.0 that includes an improved simulation of 
moderate pt junction loop with a factor Fp^—?>. 

A similar conclusion can be drawn from the results obtained for p/p ratio for central (0-10 
%) Au+Au coUisions at 200 GeV presented in Fig. 7 and for the p/Tr"*" ratio (not presented 
here). These results show the significant contributions of protons and anti- protons yields 
to the total particle composition in this moderate pt region (2 < p^ < 5 GeV/c). An 
alternative interpretation of the observed increase with centrality is provided by parton 
recombination and fragmentation model [12] while the hydrodynamics [43] and thermal 
model calculations predict that p/vr" ratio exceeds unity for central collisions. 



D. Nuclear modification factor versus pr 

In order to better quantify the particle composition at moderate pr we investigate the 
binary collision scaling of Pt spectra for charged pions and protons (anti-protons). Figure 8 
shows the predicted nuclear modification factors Raa and Rep for sum of protons and anti- 
protons (p-|-p)/2, and charged pions in comparison with available PHENIX data [23]. 

Rep is defined as the scaled yield ratio at different centrality such as the ratio of central 
to peripheral yield: 

Yield{Central) / < Ncoii{Central) > , . 

^'P^P^> = Yield{Periph.)/ < Neoii{Periph.) > 

where, Yield — {1 / N events) {^/'^'^ Pa.) {d'^N / dp ±_dy) and < Ncoii > is defined as above. 

The scaling behaviour of the pions is different from those of the sum of protons and anti- 
protons. The pions yield scaled by Ncou in central events is strongly suppressed compared to 
pp reactions {Raa) and to peripheral events {Rep)- The hadron production in HIJING/BB 
v2.0 is mainly from the fragmentation of energetic partons. Thus, the observed suppression 
in central collisions may be a signature of the energy loss of partons during their propagation 
through the hot and dense matter (possibly QGP) created in the collisions, i.e jet quenching. 

HlJlNG/BB v2.0 predicts that the sum of protons and anti-protons (p+p)/2, scale wells 
with the number of binary collisions (A^coh)- On the other hand the discrepancy seen at 
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Pt < 1-5 GeV/c indicates a sizeable contribution from radial flow. Similar trend are observed 
in A, Kg and K"^ measurements by the STAR Collaboration [65]. It has been recently shown 
that the competition between recombination and parton fragmentation at moderate pt may 
also explain [12] the observed features. 

IV. SUMMARY AND CONCLUSIONS 

The failure of the previous implementation of baryon junction loops in HIJING/BB 
vl.lO to reproduce the observed pt enhancement of anti-baryons and baryons motivated us 
to construct a new version, HIJING/BB v2.0. 

The new physics ingredients implemented in HIJING/BB v2.0, concentrate on modifica- 
tions of the hypothesized junction JJ loops algorithm as well as adding a phenomenological 
simulation of collective transverse flow. These modifications make a significant improvement 
in the full event Monte Carlo description of a large set of observables for p and p. The new 
version can account now for many features of the baryon anomaly region at moderate pr, 
as well as for characteristic stopping observables at ^/sNN— 200 GeV in Au+Au coUisions. 

A simultaneously absence of suppression for baryons up to pr=4-5 GeV/c and the en- 
hancement of the p/tt ratios at moderate px which is a challenge for many theoretical frame- 
work is well described within HIJING/BB v2.0 with shadowing and jet quenching effects, 
for an energy loss parameter dE/dx =1 GeV/fm (for quark jet) and a constant phenomeno- 
logical factor Fpj,—3. One of the remaining discrepancies is the energy loss per participant 
nucleon and baryon rapidity measurements at forward rapidity ( y > 3), which will require 
futher analysis. 

While this new version HIJING/BB v2.0 gives a good description of a large body of 
data it still cannot reproduce the transverse mass spectra of kaons for which the integrated 
yield is well predicted, but the model has no mechanism to account for kaon radial flow. In 
string fragmentation phenomenology, the strong enhancement of strange particle observables 
require strong color field effects (SCF) [58,59]. The full understanding of the production of 



13 



strange particles in relativistic heavy-ion collisions [66] remains an exciting open question. 
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FIG. 1. (Color online) Comparison of HIJING vl.37 (left part) and HIJING/BB v2.0 (right 
part) predictions for the nuclear modification factor {Raa) for central (0-10%) Au+Au collisions 
at y/sNN=200 GeV. The results are with (solid histograms-yqs) or without (dashed histograms-nqs) 
shadowing and quenching effects included. The label f3 stands for model calculations assuming 
Fpj, = 3. The data are from PHENIX [62]. Only statistical error bars are shown. The error bars 
at pt=2 GeV and pt=3 GeV, include systematic uncertainties. Equivalent error bars on the other 
points have been omitted for clarity. 
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FIG. 2. (Color online) HIJING/BB v2.0 (upper panel) and HUNG vl.37 (lower panel) pre- 
dictions of transverse mass distributions for positive pions (vr+) and protons (p) (left part), and 
negative pions (vr^) and anti-protons (p) (right part). The solid and dashed histograms have the 
same meaning as in Fig. 1. The data are from PHENIX [22]. The error bars show statistical errors 
only. 
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FIG. 3. (Color online) Calculated rapidity dependence of a) protons, b) anti-protons, c) 
net-protons, and d) anti-proton to proton ratio. The solid and dashed histograms have the same 
meaning as in Fig. 1. The dotted histograms correspond to RQMD v2.4 model predictions. The 
data, corrected for weak decays, are from PHENIX [23] and BRAHMS [28] . The errors bars shown 
includes both statistical and systematic uncertainties. 
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FIG. 4. (Color online) Comparison of HIJING/BB vl.lO (dashed lines) and HIJNG/BB v2.0 
(soild lines) predictions for the centrality dependence of proton (part a) and anti-proton (part b) 
yields at mid-rapidity. The data are from PHENIX [23] . The errors shown includes both statistical 
and systematic micertainties. 



23 



120 
I 80 

o 

^ 60 
% 40 
20 





(a) Au+Au, 200A GeV (cm.) 

Central (0-5%) 

• PHENIXdata 

A BRAHMS dataj 
<(5yeB> = 1 .75 
<(5y,3,> = 1 .67 
<(5y,,> = 2.57 



(b) Au+Au, 200A GeV (cm. 
Central (0-5%) 



) 



Ebb=42.6 (GeV) 
E,3,=40.5 (GeV) 
E,,=23.0 (GeV) 




4 6 

rapidity, y 



rapidity, y 



2.5 

C 

o 
o 

I 1-5 
-*-> 



^ 0.5 



(c) Au+Au, 200AGeV (cm) 
Central (0-5%) 

BB2.0 (yqs,f3) 

HIJ137 (yqs) 

RQMD v2.4 



rapidity, y 

FIG. 5. (Color online) Model predictions for: a) the net-baryon distribution and the average 
rapidity loss < by >\\>) the total energy per net-baryon after the collisions; and c) the transverse 
energy per net-baryon for central (0-5%) Au+Au collisions at y^J/v7V=200 GeV. The solid and 
dashed histograms are the results obtained within HIJING/BB v2.0 and HIJING vl.37, respec- 
tively. The dotted histograms are the RQMD v2.4 predictions. The data are from BRAHMS [28]. 
The errors bars include both statistical and systematic uncertainties. 
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FIG. 6. (Color online) Model predictions for p/vr" ratio in central (0-10 %) Au+Au collisions 
at 200A GeV (upper part) and peripheral (60-90 %) Au+Au collisions (lower part). The solid and 
dashed histograms have the same meaning as in Fig. 1. In figures b and d, the dotted histograms 
are the predictions of HIJING/BB v2.0 with Fy^=\ (label fl). The data are from PHENIX [23]. 
The error bars include systematic uncertainties. 
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FIG. 7. (Color online) Comparison of HIJING vl.37 (left) and HIJING/BB v2.0(right) model 
predictions for p/p ratio versus pT for central (0-10 %) Au+Au collisions at 200A GeV. The solid 
and dashed histograms have the same meaning as in Fig. 1. The data are from PHENIX [23]. The 
error bars include systematic uncertainties. 
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FIG. 8. (Color online) HIJING/BB v2.0 predictions for binary-collision scaled nuclear modifi- 
cation factor Raa (part a) and Rep (part b) for (p-|-p)/2 (solid histograms) and charged pions 
(dashed histograms) in central (0-10 %) Au-|-Au collisions. The data are from PHENIX [23]. The 
error bars are statistical only. 
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